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This paper develops energy-transfer modelling of active thermoacoustic engines resonantly controlled on 
boundary for amplification of power rating toward satisfaction of renewable industry. Therein the wave 
equation of thermoacoustic dynamics in resonators with non-uniform media and boundary actuations is 
derived and then turned into a least-action principle. With this least-action principle, we obtain the gov¬ 
erning equation of longitudinal resonators with spatially variant cross-section areas to investigate how to 
shape the resonator for boosting piston stroke and power-transmission efficiency. It is followed by spa- 
tiotemporal transfer-function modelling that functionally represents the dynamics and interprets the 
boundary actuations into internal inputs. This helps formulate the overall dynamics into feedback-inter- 
connection between the thermoacoustic dynamics in the resonator and the mechatronic dynamics of the 
alternative current generator, so that synthesis of feedback systems can be applied to design the entire 
engine. Transfer-function modelling following least-action principle leads to the conservation law of 
thermoacoustic storage, which figures out engine cycles, the most fundamental principle in designing 
active thermoacoustic engines. Based on such feedback realization, digital signal processing is pro¬ 
grammed to numerically assess power ratings of active designs. 

© 2014 Elsevier Ltd. All rights reserved. 


1. Introduction 

Of facilities generating electricity from sunlight, the 
thermoacoustic engine is of much less ratio of manufacturing cost 
to power rating than photovoltaic panels [1,2], solar-thermal 
devices, and other heat engines. Its extremely simple mechanism 
enables us to build solar power plants in poor lands and home roof¬ 
tops at affordable prices. Besides, thermoacoustic engines under 
resonant control can have power efficiency even better than Carnot 
cycle efficiency that merely exists in an ideal heat engine. Inas¬ 
much as passive thermoacoustic engines, for examples in [3-10], 
poorly generate electricity, this work aims to investigate active 
thermoacoustic engines with power ratings toward the require¬ 
ment of renewable energy industry. 

The active thermoacoustic engine under design is assembled by 
alternative current generator, Rijke tube, and radiation-to-heatflux 
actuator, as shown in Fig. 1. The radiation-to-heatflux actuator is 
comprised of crystal liquid photo valve and ultra-conducting por¬ 
ous media adjacent to each other. It consumes very little power 
to actively switch the heat-flux source at resonant frequencies 
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for amplification of power rating. At the head end of the Rijke tube, 
the heat-flux resulting from sunlight gives rise to entropy rate that 
changes the fluid density there, which in turns initiates acoustic 
velocity. That is, the heat-flux excitation expanding and contract¬ 
ing the fluid acts as a source of acoustic wave that rides on 
Rijke-tube media and arrives at the load end to push and pull the 
piston of alternative current generator. 

A frequency generator embedded therein is programmed to 
online keep the switching frequency of the buck chopper that drives 
the photo valve on one of resonant frequencies, which prompts 
standing waves up to resonance and thus increases mechanical 
power enormously. Besides, the Rijke tube is longitudinally shaped 
to make power transmitted from the heatflux actuator to the cur¬ 
rent generator be with high pressure and low velocity. Such high- 
pressure AC power transmission reduces the viscosity effect inside 
the Rijke tube and increases the piston stroke, analogously to adjust 
the gear ratio in mechanical transmission. The stackers inside the 
Rijke tube are for internal heating of working gas to maintain 
higher-frequency resonance for larger power rating. For the piston 
at the tail end to move back and forth, some clearance between the 
Rijke tube and the piston surface produces acoustic leakage and 
friction on the right boundary. We model this practical situation 
by Robin boundary condition. 
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Fig. 1. Schematic design of resonant thermoacoustic engines. 


The thermoacoustic dynamics inside Rijke tubes has been inten¬ 
sively studied in recent years, for examples, based on impedance 
analysis [6-10], computational fluid dynamics [11-14] or spatio- 
temporal transfer-function [15]. In fact, these studies mostly fall 
heirs to the thermoacoustic dynamics in the field of combustion 
instabilities, for examples in [16-20]. No matter in thermoacoustic 
engines or in combustion instabilities, the kernelled dynamics is 
concerned with the process how heat perturbation excites acoustic 
motions, known as Rayleigh Criterion [15,21-24]. In this paper, the 
Rayleigh criterion is further formulated into a principle of least 
action, named here by Lagrangian thermoacoustic dynamics. It is then 
employed to derive the wave equation governing energy transfer 
through Rijke tubes with spatially variant cross-section areas. 

With Lagrangian formulation, cylinder’s pressure-velocity 
amplitudes are thus capable of being longitudinally adjusted mainly 
for energy efficiency in power transmission. As the cross-section is 
nearly exponentially narrowed along the longitudinal direction, 
power is transmitted from the head end with larger pressure and 
tiny velocity to the load end with smaller pressure and required 
velocity. Low-velocity AC-power transmission reduces the viscosity 
effect inside the Rijke tube and thus increases energy efficiency. 

The wave equation with non-uniform media in conjunction with 
inhomogeneous boundary conditions is then transformed into a 
single spatiotemporal transfer function, which is an expanded 
application of the functional tool known as nD transfer function 
models [25-29]. Therein the inhomogeneous boundary conditions 
arise from the von Neumann source of heat-flux at the head end 
and the Robin source of piston-speed at the tail end, both of which 
are transformed into Dirac internal sources by way of Laplace- 
Galerkin transform. This renders the capability of representing 
the overall dynamics into feedback-interconnection between the 
dynamics of the Rijke tube dynamics and the alternative current 
generator, and thus of transforming the engine design into synthe¬ 
sis of feedback systems. Moreover, digital signal processing can be 
programmed based on this feedback realization to offline or online 
assess the power rating. 

A new conservation law is inferred from transfer-function trans¬ 
formation of Lagrangian thermoacoustic dynamics, named by con¬ 
servation law of thermoacoustic storage. It figures out engine cycles 
and thus is the most fundamental principle in designing active ther¬ 
moacoustic engines. For example, flow leakage through the gap 
between piston surface and chamber interior happens in real 
designs, which results in a Robin source to the load end. Fortu¬ 
nately, we can estimate the level of flow leakage with the conserva¬ 
tion law of thermoacoustic storage through identification of Robin 
coefficients therefrom. As we know, Robin coefficients can only be 
determined through system identification [30,31], and this paper 
provides a method to identify the Robin coefficients of flow leakage. 


2. Lagrangian thermoacoustic dynamics 


Let Q be a Rijke tube containing non-uniform gas medium with 
time-invariant, spatial variant entropy s, pressure p, temperature 
T, and density p. Imposed on the medium is the thermoacoustic 
motion with entropy fluctuation s, acoustic pressure p, and acous¬ 
tic velocity u. The following assumptions are made: 

(1) The medium is an ideal gas. 

(2) The heat transfer and viscosity in the field of fluctuation is 
neglected. 

(3) The medium is stationary, i.e., ii = 0. 

The dynamics of thermoacoustic motion was usually derived 
from the first-order perturbation of fluid conservation equations 
[32 . Application of thermoacoustic state transformations in [33] 
to this set of perturbation equations yields a new phrase as follows: 


ds 

dt 


-Vs • u + 


d 

PT’ 


du 

dt 



( 1 ) 

( 2 ) 


= -kV ■ (pu) + (7 - 1)(Q - pTVs ■ u), k = yRT , (3) 

where Q stands for the source of heat-perturbation rate. Therein the 
parameters R, y and k are the gas constant, the specific heat ratio 
and the compressibility of the gas medium, respectively. 

Eq. (1) implies that spatial variation of mean entropy s on the 
border of any differential volume moving with bulk velocity 
u(u = 0) produces entropy flux -Vs ■ u through the border. The 
total entropy-rate of the differential volume is the entropy flux 
-Vs ■ u plus the entropy change Q/pT provided by the heat-per¬ 
turbation rate Q. Eq. (2) is directly from the Newton's second 
law. Eq. (3) is based on the state transformation: 


P_ = P + 1_ 

yp P <V 


( 4 ) 


Combined by the law of continuity: 

f=-V-(pu) (5) 

Eq. (4) results from the linearized Gibbs equation in conjunction 
with the linearized state equation for an ideal gas. It implies that 
“pressure change can be produced by fluid compression or entropy 
fluctuation" [33]. Substituting Eq. (1) and the temporal differential 
of Eq. (5) into Eq. (4) yields Eq. (3). 

Then, substituting Eq. (2) into the temporal differential of 
Eq. (3) yields 

cPp __ 2 , du . DQ_ 

-^- K V p=-(y-l)pTVs. —+ (y-l)— (6) 

Defining the thermoacoustic momentum ijj by the temporal inte¬ 
gration of acoustic pressure p, that is, 


p = ^; pii=-Vi/q (7) 

We rephrase Eq. (6) by 

S-(KV + (y-l)TVs)-V* = (y-l)Q, in O. (8) 

dt 

Firstly, Eq. (8) can be transformed to a wave equation free of 
damping: 
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d 2 i/> 

W 


_J_ V .(e(W P V^) = ( 7 -l)Q, 


(9) 


where the s 0 is the mean-entropy at an arbitrarily chosen location 
x 0 inside or outside the chamber Q. Therein the parameter C p is 
the constant-pressure heat capacitance. The derivation from Eq. 
(8) of Eq. (9) is as follows. Eq. (8) is identical to 


g> 2 i/> 

W 


- K[ Vl/f + 


Vs 

Cp 



(y-m, ^ 


yR 

7-1' 


In juxtaposition of this equation with the canonical form of 
wave equation, such as Eq. (9), it is found that the existence of 
an integration factor a such that 

er((Vs/C p ) • Vi// + V 2 <//)) = V • (erVi/O 


is necessary to transform Eq. (8) to Eq. (9). The above equation 
implies Vo/a = Vs/C p , that is, o = exp f c Vs ■ dr), where C is 
an arbitrary path for line integral, wherein the line integral is inde¬ 
pendent of path since the mean-entropy s belongs to thermody¬ 
namic state. Therefore such a factor o exists, which equals 
a = e (s ~ s °>/ c p. With this found factor, it can be verified that Eq. (8) 
is identical to Eq. (9). 

Secondly, the state equation combined by the first law of ther¬ 
modynamics [33] equates 


e (i-s 0 )/Cp _J_( _P_\ 

To \PoJ 


( 10 ) 


where f 0 and p 0 are the mean temperature and pressure, respec¬ 
tively, at the same reference location x 0 . Substitution of Eq. (9) for 
Eq. (10) yields 


that calculus of variation upon Eq. (14) also leads back to the wave 
equation of Eq. (12), named here by Lagrangian thermoacoustic 
dynamics. 

3. Spatiotemporal transfer function modelling 

Fig. 1 shows the thermoacoustic engine schematics under design 
that comprises that comprises the Rijke chamber, the radiation-to- 
heatflux transducer, the photo valve, and the alternative current 
generator. Here we consider the thermoacoustic dynamics to be 
longitudinal in the Rijke tube with spatially variant temperature, 
pressure and cross-section area. Moreover, the flow leakage and 
boundary friction are also concerned. 

The heat-flux excitation q (W/m 2 ) on the head end (x = 0“) acti¬ 
vating thermoacoustic process inside the Rijke chamber is equiva¬ 
lent to a pointed heat generation q< 5(x) between x = 0“ and x = 0 + 
with heat isolation at x = CT, based on the localization of energy 
conservation. Here d is the Dirac delta distribution. With Eq. (13), 
the Rijke chamber stores Rayleigh’s Kinetics K R and Rayleigh’s 
Potential P R being, respectively, 

r i r i 

2K r = / R a i// 2 Adx and 2 P R = / KR a x//' 2 Adx , (15) 

Jo Jo 

where A stands for the cross-section area of the Rijke chamber that 
is spatially variant. The Lagrange formulation of the Least-action 
principle of Eq. (14) upon Eq. (15) yields the following longitudinal 
dynamics: 


d V 
W 




(y-m 


( 11 ) 


It can be verified that f 0 and p 0 in Eq. (11) can be removed 
without any change, which is in line with the fact that the dynam¬ 
ics has to be independent of the choice of reference location x 0 . 
Defining the dimensionless Rayleigh factor by R a = (p/p 0 ) (1 
we express the thermoacoustic dynamics by 


Tjpr~ jjrV ■ {KRaVi//) = (y - 1)Q, (12) 

where k = yRT denotes the compressibility. Without losing gener¬ 
ality, the reference pressure p 0 can be chosen as the environmental 
pressure. 

To realize the dynamic behaviour described in Eq. (12), Rayleigh 
[34] claimed the analogy of acoustic pressure p to mechanical 
velocity and heat-rate source Q to excitation force, which has been 
largely investigated in combustion instabilities earlier or thermoa¬ 
coustic engines later as the Rayleigh criterion. To quantify this 
analogy based on Eq. (12), we define two internal scalars: Ray¬ 
leigh’s Kinetics l( R , Rayleigh's Potential P R and one line integral Ray¬ 
leigh’s Work W R respectively by 


2 K r = j^R a 00 dV, 2P R = £ KR a \Vxl/\ 2 dV, 

dW R = f dxf ■ R a (y - \)Q_ dV. 

Jn 


(13) 


With calculus of variation, Eq. (12) leads to a least-action 
principle: 

0= [ 2 5I( R ~SP R + dW R dt, (14) 

Jt, 


with respect to the variation of Rayleigh’s displacement 5i//, where <5 
represents the operator of functional variation. It can be verified 


in addition to piston vibration at x = K Since the reference pressure 
po has been chosen as the environmental pressure, the Rayleigh fac¬ 
tors on both boundaries are of unity, i.e., R a (0) = R a (H) = 1. Reduction 
of non-uniform medium to be spatially invariant and replacement 
of heat-flux excitation by acoustic driver, Eq. (16) appears to be 
Webster’s horn equation which is well-known in the acoustics of 
loudspeakers and flared wind and brass instruments. 

Eq. (16) facilitates the investigation of the thermoacoustic pro¬ 
cess happened on the head end between x = 0~ and x = 0 + . The 
hyperbolic Eq. (16) has signified that the acoustic pressure p in 
Rijke chamber is spatially continuous at x = 0, i.e., there is no pres¬ 
sure jump from x = (T to x = 0 + . Thus, the integration of Eq. (16) 
between x = CT and x = 0 + yields the left boundary condition to be 

CpTi/i' = -q, at x = 0 + , (17) 

where C p = Ry/(y - 1). Based on the definition in Eq. (7), Eq. (17) 
can be rephrased by 

u = Cll S- atx = 0 + (18) 

7 Po 

Eq. (18) reveals that the heat-flux q initiates an acoustic velocity 
on the left boundary, but keeps the acoustic pressure continuous 
on the boundary. 

The dynamic process on the head end as indicated by Eq. (18) 
can be realized by tracing back to the Chu’s formula of Eq. (4) 
and the continuity law of Eq. (5) as follows. In between x = CT 
and x = 0 + the heat-flux excitation q gives rise to entropy rate s 
that changes the fluid density, as indicated by Eq. (4), which in 
turns initiates acoustic velocity, as indicated in Eq. (5). That is, 
the heat-flux excitation expanding and contracting in fast-time 
scale the fluid on the boundary x = 0 acts as a source of acoustic 
wave. The acoustic wave rides on Rijke chamber media and arrives 
at the load end at x = t to push and pull the piston, generating 
mechanical work. 
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As for the boundary condition at the load end, it is nominally 
taken to be 


i y = -pv, at x = £ , 

with the assumption that the acoustic pressure is continuous at 
x = t and there is no heat-flux excitation at x = tt, where v stands 
for the piston speed. If we take flow leakage and boundary friction 
into consideration, the right boundary condition belongs empiri¬ 
cally to the Robin type: 

Gap + P4>' =-pv, atx = t~, (19) 

where a > 0, p > 0. Here the leakage factor p concerns the ratio of 
closed area to open area on the load-end boundary, which is to be 
identified in real design. Moreover, the piston dynamics chosen 
for representation is governed by 

mv-\-bv = \pt, (20) 


wherein m stands for the inertia of piston and b comprises the 
boundary and load damping. For succinct notation, i //(x, t ) may 
appear as l j/ x (t) somewhere in this paper and similarly to other 
variables. 

The collection of Eq. (16) to Eq. (20) describes the thermoacou¬ 
stic engine dynamics as 


W 


1 d 
R a A dx 



0, in (0, £)\ 


( 21 ) 


to describe f s modal decomposition in the s-domain. Accordingly, 
the inverse Laplace-Galerkin transform is the composite of the 
inverses of those two ID transforms, i.e., MT 1 = = if -1 

'f, '. Explicitly, 

f{x,t) = jtr'lFds)] = 

With the above theoretical background, in the next paragraph 
we will perform the Laplace-Galerkin transform on Eqs. (21)-(24) 
to derive the 2D transfer-function modelling of the engine 
dynamics. 

At the load end x — l, substitution of i.fx + Pf\ = 0 for Eq. (23) 
x( <i>i) yields 

(M - = (4>x/P)pv at x = i, (28) 

and the integration by parts yields 



+ KAi/fa I x= o + ~ M') \x=e- 

Based on Eq. (22) and (28), the above equation becomes 

J‘-JL^KRaA^ ■ <t>,dx = - (y - 1 )qA 0 ■ UO) 

+ pr'yPoAiV ■ <p x (e), (29) 


$>(*) [ e«F(X,s)ds. (27) 

icA J r 


CpTif/' = -q, at x = 0, 

oap + Pf' = -pv, at x = £(a > 0, p > 0); and 


( 22 ) 

(23) 


where p 0 denotes the environmental pressure, and therefore R a ( 0) 

= RaW = 1. 

With Eq. (29), performing Laplace-Galerkin transform on Eq. 
(21) yields 


mv+ bv = \ji, at x = l, 


(24) (s 2 +2)f'(2,s| = (y- 1)A 0 ^(0) ■ q(s)-p 0 A,p 'y^iW ■ Hs), (30) 


wherein the Rijke-tube dynamics of Eqs. (21)-(23) is feedback-cou¬ 
pled with the represented load dynamics in Eq. (24). 

The Rijke-chamber dynamics in Eqs. (21)-(23) involves the 
elastic stiffness operator : 


jfV/; = - 


1 d 

R„A dx 



the domain of which is 


D(jT) ee {</> : = 0 at x = 0; etc/) + P4>' = 0 atx = £}, 

where a > 0, and p > 0. The stiffness operator ,jT is self-adjoint and 
positive-definite under the inner-product: 

((p,cl>)= f R a (x)A(x) <p*(x)<p(x)dx. (25) 

Jo 

Denote by A the eigenvalues set of the stiffness operator jP\ 
corresponding to the eigenfunctions set <h = {0i} ieA , i.e., 

.if r/); = /</). in [0, £], 4>; e D(jT) for V2 e A. 


where W(J,s) = if[4i(x, t)].q(s) = iP[q(t)], and v(s) = £f[v(t)]. The 

2D transfer function of the Rijke-chamber dynamics of Eqs. (21)- 
(23) is thus obtained to be 


= (y - 1)A 0 


<k(Q)q(s) 

s 2 + X 


PoA,P l y 


4>i(e)Hs) 

s 2 +1 


(31) 


feedback-coupled with the load dynamics: 

= -b^’ (32) 


which is obtained by performing Laplace transform on Eq. (24). 

Eqs. 31,32 constitute the spatiotemporal transfer-function mod¬ 
elling of the thermoacoustic engine. This represents the engine 
dynamics as a feedback-interconnection of the chamber dynamics 
in Eq. (31) and the load dynamics in Eq. (32), which transforms 
the overall design into synthesis of feedback systems. Moreover, 
the feedback-coupled transfer-functions can be realized in the state 
space for easily developing digital signal processing for online or 
offline computation of engine performance. 


Accompanied by j^’s self-adjointness is that its eigenfunctions 
set is real and orthonormal, and followed by :/f's positive-defi¬ 
niteness is that all of its eigenvalues are positive, X > 0 for 
Vi e A. Moreover, the countable eigenfunctions set is a complete 
basis of L 2 ([O.t']), which follows that the inverse of .if is a bounded 
operator and closed in the Banach space L 2 [0.£]. 

With respect to the eigenfunctions set <£ and the inner-product 
of Eq. (25), a 2D transform .if, named by Laplace-Galerkin trans¬ 
form, is then defined by the composite of Laplace transform and 
Galerkin transform Zf = i.e., 

/*oo nX=t 

F{X,s) = if\f(x,t)]= / e- st R a (x)A(x)<l> l (x)f(x,t)dxdt, (26) 

Jt =0 Jx =0 


4. Conservation of thermoacoustic storage 

Performing the inverse Laplace-Galerkin transform ' on the 
spatiotemporal transfer-function of Eq. (31) yields 

w it S ( KM s) - 

-Pop-'y v{t)Ai52<l> x (£)<l> x (x), (33) 

AeA 

with homogeneous boundary conditions: i//(0) = 0 and ai/'(^)+ 
ptpfe) = 0. Eq. (33) equals 
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w~^A^ {KRaA ^ = (y- www -p°p x y v(t)5(x -i). 

(34) 


Eventually, the wave equation of Rijke tube in Eqs. (21)-(23) is 
mathematically identical to Eq. (34) with homogeneous boundary 
conditions. In this way, spatiotemporal transfer-function enables 
us to realize inhomogeneous boundary conditions as internal 
sources in conjunction with boundary homogeneity. 

The operator form of Eq. (34) becomes 

{J/ij/ + JTt/0 (x, t) = (y - l)q(t)<5(x) - Po/T 1 ? v{t)5(x - £), (35) 


where the mass operator M is the identity operator, and the stiff¬ 
ness operator :ZZ that has been detailed in Section 3 is 


yff = - 


1 d 

R a A dx 



(36) 


with the domain of operation 

Dpf) = {i/r: <t>'(0) = 0, = 0}. 

Both operators are positive-definite and self-adjoint under the 
inner product (■,•) defined in Eq. (25). 

Define a new quantity E, named by thermoacoustic storage, to be 

2 E = (i/r, Ji^i) + (i If, JT i \i). (37) 

The thermoacoustic storage E is positive-definite, since the 
operators M and .Zf are positive-definite under the inner-product 
defined in Eq. (25). It becomes zero if and only if \j/ = 0 and \]/ = 0 
almost everywhere. The operators M and r/f are self-adjoint, so 

ft = j Jl ^) + ('/'> J/ ^) + W'I') + (•A, W¥i) 

= (l ]l, + (l/q JTl/t) = (l/q + C/C\j/). (38) 


Based on Eq. (35 ), the above equation becomes 

E = Ao(y-\)qp 0 -p 0 A t r'yvp„ (39) 

where p = \j/ denotes the acoustic pressure. 

Further to define E T by scaling the thermoacoustic storage E to 
be 


E t = -4- (40) 

ypo 

so that Eq. (38) becomes 

E T =A 0 { -^Mp 0 -r%v Pl ( 41 ) 

yp o 

Based on Eq. (18), the first term of the right side of Eq. (40) 
equals A 0 u 0 p 0 , where u 0 is the acoustic velocity at the head end 
x = 0. Then Eq. (40) can be written by 

E t = AqUqPo - p-'AtUi p e , 

which reveals that E T defines the mechanical energy stored in the 
thermoacoustic dynamics. The mechanical energy E r extends the 
mechanical part of thermoacoustic energy for non-uniform pressure 
and temperature from that in uniform mean-flow conditions that 
has been defined in literature of thermoacoustic dynamics, for 
example, in [33], 

Suppose the boundary heat-flux q entering into the chamber is 
of resonant frequency co, in steady state the thermoacoustic stor¬ 
age E will be conserved over a period T = 2nl<±>, that is, 

A t [ v(t)Pi(t)dt = Mo (y - ^ [ q(t)p 0 (t)dt, (42) 

Jt yp o Jt 


which is to be identified by Eq. (41 ). Furthermore, the integrals of 
Eq. (41) can be written as the integration over a cycle of thermoa¬ 
coustic momentum i jz, that is, 

[ Vd^ = p^ t 1 ! [ qd^ 0 , (43) 

Jw >\i yp o Js 

where W and S denote the work cycle and source cycle, respectively. 

Over a cycle, the points of (i/q, v) depicts a closed region, the 
area of which indicates the mechanical energy generated or lost 
over a circle. If the points of (i/q, v) direct clockwise along the work 
cycle, it means that positive power is gained. Similarly, if the points 
of (\j/ 0 . u 0 ) direct clockwise along the source circle, the useful power 
excited by the boundaiy heat-flux is positive. As such, Eq. (42) 
defines the engine cycles for active thermoacoustic engines, as 
shown in Fig. 2, where the upper parts of the regions enveloped 
by (i j/,, v) and {\j/ 0 , u 0 ) indicate the expansion stroke, and the lower 
parts indicates contraction stroke. Unlike Stirling engines that have 
four strokes including heating and cooling strokes, this thermoa¬ 
coustic engine has only two strokes so that they can be fulfilled 
by a single cylinder. 

Engine cycles in Eq. (42) reveals the following design principles: 

(1) The heat-flux excitation can be on-line controlled to be of a 
frequency identical to one of the resonant frequencies to 
arrive at the maximum power generation. 

(2) The power rating of engine is proportional to the square of 
(y - 1 )ly, concerned with the working gas. 


xlO 3 



Ay/ {kg ■m/sec ) 







Source cycle 
Work cycle 
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where the left side is the mechanical WOlk geneiated ovei a peiiod. pjg 2. Comparison of engine cycles at resonance with those without resonant 
There, the leakage factor fl< 1 quantifies the loss of work in practice, control. 
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(3) The ratio of the heat-end area A 0 to the load-end area A f 
plays like a gear ratio in power transmission adaptable to 
the mechatronic load. 

(4) The engine is getting more powerful along with higher 
working frequency of resonance that is usually constrained 
by control rig and working temperature. 

5. Performance review 

Let the thermoacoustic engine of Fig. 1 be with solar heat-flux 
excitation q = 500 W/m 2 , length of chamber f = 3m, cross-section 
areas of head end and load end A 0 = 7 m 2 and A t = 0.7 m 2 , respec¬ 
tively, the working gas of specific heat ratio y = 1.4 and universal 
constant R = 287.058 J/kg K, leakage factor p = 0.9, and working 
temperature T = 300 K in average. The mechatronic load is repre¬ 
sented by the inertia of 10 kg series-connected by the damping 
of 4 nt s/m. 

Fig. 2 shows working and source cycles coordinated by Ai/c and 
u, wherein the first plot is without resonant control and the second 
plot is at resonance of the first-mode frequency /] = 14.64 Hz. 
Therein both cycles are directed clockwise, meaning that the 
engine is doing positive work on the load end and the boundary 
heat-flux is exciting positive useful energy. The closed region 
enveloped by the work cycle has the area A w meaning the mechan¬ 
ical work done by the engine over a cycle. Similarly, the enclosed 
area by the source cycle A s is the useful energy trigged by the 
heat-flux q. The first plot verifies the conservation law presented 
in Eq. (42): 

Aw = flAs , 

where p = 0.9 is the leakage factor indicating the energy loss in 
practical design. We capture this energy-loss effect by modelling 
Robin-type boundary condition upon the load end. In both cycles, 
the I—II phases indicate the expansion strokes, and the III—IV phases 
indicate the contraction strokes. 

As indicated by the second plot, while the working frequency is 
approaching to resonant frequency, the engine cycles are getting 
extreme in shapes and thereby in power rating, which demonstrate 
the main difference between active design and passive design. For 
clearer view of the engine cycles at resonance, Fig. 3 shows its 
responses in time domain, where p 0 , p e stand for the acoustic pres¬ 
sures at the head end and at the tail end, respectively, and v for the 
piston speed. With the first-mode resonance combined by the ratio 
A,/A 0 = 1/10, the power is transferred from the head end with 
A 0 ij/ 0 = 1837.21 kg m/s and u 0 = 0.0014 m/s of peak values to the 
load end with Ap\i t = 67.80 kg m/s and v= 6.78 m/s of peak values. 



Fig. 4. Piston strokes at different resonant frequencies and head-to-tail area ratios. 

Analogously, AC power is transmitted by much larger pressure and 
tiny velocity to the load end with proper pressure and velocity, so 
that the viscosity effect inside the Rijke tube is significantly 
reduced. That is, the ratio of the heat area A 0 to the load area A, 
functions like AC-power transformer in power transmission. Mean¬ 
while, it results in long piston-stroke, which makes possible in 
mechanism design large gear-ratio transmission from mechanical 
thrust to electricity generation, and thus increases the power-gen¬ 
eration efficiency. 

Figs. 4 and 5 record the responses of piston stroke and power rat¬ 
ing, respectively, as the photo valve is switched at the zero, first, 
second and third -mode frequencies and the area-ratio A 0 /A, is 
set from 1 to 10. Fig. 4 shows that the piston stroke is boosted by 
increasing the area-ratio A 0 /A,. For example, at first-mode reso¬ 
nance, the piston stroke d s is boosted from 0.12 cm to 16.75 cm as 
the area-ratio A 0 /A t is changed from 1 to 10. Due to the boundary 
friction arising from flow leakage, resonance at higher frequencies 
yields shorter piston-stokes with a specified area-ratio. As a result, 
the area-ratio in thermoacoustic engine plays analogously as gear- 
ratio in power-train transmission to regulate force-velocity. 

Fig. 5 contains the Bode plot of Po/Patnv which shows the ampli¬ 
tude of head pressure p 0 and its phase difference 0 between the 
head velocity u 0 , across the span of working frequencies. The effect 
of resonance on power rating can be easily examined by this Bode 
plot. In the magnitude plot, four resonant frequencies are marked: 
/o = 5.37 Hz,/] = 14.64 Hz,/ 2 = 24.30 Hz and/ 3 = 34.29 Hz, at each of 



(/n/sec) 



t (sec) 



frequency ( rad /sec) 


Fig. 3. Steady-state responses of head pressure, tail pressure and piston velocity at Fig. 5. Frequency response of the head pressure and power rating assessment at 
resonance. different resonant frequencies. 
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which the real power POW generated by sunlight is marked above 
each peaks. The real power POW is calculated by 

POW = pU 0 Po cos 0, 

where U 0 and P 0 are root-mean-square values of u 0 and p 0 , respec¬ 
tively, and the phase difference 0 becomes the power factor in AC 
power. Generally, the power rating is getting larger along with res¬ 
onance at higher frequencies, but the piston stroke is getting 
shorter. 

Suppose that the head area is A 0 = 7 m 2 and the efficiency of 
electricity generation is 60%, resonant control working at the sec¬ 
ond-mode frequency can harvest 102.07 W. A suburban family 
needs almost 1.0 KW to constant power it up. For any homes pos¬ 
sessing more than 70 m 2 rooftop or land spaces, the resonant ther¬ 
moacoustic engine under investigation is the best candidate served 
for solar energy generation to power them up or turns electricity 
metres backward. 

6. Conclusion 

This paper develops the energy-transfer modelling of active 
thermoacoustic engines and performs assessment of power rat¬ 
ings. We find that the power rating and piston stoke can be 
remarkably amplified through resonant control of boundary heat- 
flux and by shaping the resonator, respectively. Such findings are 
by way of Lagrange thermoacoustic dynamics, spatiotemporal 
transfer-function, and the conservation law of thermoacoustic 
storage. 

Acknowledgement 

National Science Council of Taiwan financially supported this 
work for three years. The authors here express their utmost thanks 
to NSC Control Group. 

References 

[1] Shen C, He Y, Li Y, Ke H, Zhang D, Liu Y. Performance of solar powered 
thermoacoustic engine at different tilted angles. Appl Therm Eng 
2009;29:2745-56. 

[2] Wu Z, Dai W, Man M, Luo E. A solar-powered traveling-wave thermoacoustic 
electricity generator. Sol Energy 2012;86:2376-82. 

[3] Backhaus S, Swift GW. A thermoacoustic-Stirling heat engine. Nature 
1999;399:335-8. 

[4] Garrett S, Backhaus S. The power of sound. Am Sci 2000;88:516. 

[5] Swift GW. Thermoacoustic energy conversion. In: Rossing T, editor. Handbook 
of Acoustics. Springer; 2007. 

[6] Yu Z, Jaworski AJ. Impact of acoustic impedance and flow resistance on the 
power output capacity of the regenerators in travelling-wave thermoacoustic 
engines. Energy Convers Manage 2010;51:350-9. 


[7] Yu Z, Jaworski AJ, Backhaus S. Travelling-wave thermoacoustic electricity 
generator using an ultra-compliant alternator for utilization of low-grade 
thermal energy. Appl Energy 2012;99:135-45. 

[8] Symko OG, Abdel-Rahman E, Kwon YS, Emmi M, Behunin R. Design and 
development of high-frequency thermoacoustic engines for thermal 
management in microelectronics. Microelectron J 2004;35:185-91. 

[9] Jensen C, Raspet R. Thermoacoustic power conversion using a piezoelectric 
transducer. J Acoust Soc Am 2010;128:90-103. 

[10] Smoker J, Nouh M, Aldraihem O, Baz A. Energy harvesting from a standing 
wave thermoacoustic-piezoelectric resonator. J Appl Phys 2012:111:104901. 

[11] Jung S, Matveev KI. Study of a small-scale standing-wave thermoacoustic 
engine. Proc Inst Mech Eng Part C-J Eng Mech Eng Sci 2010;224:133-41. 

[12] Matveev KI, Jung S. Modeling of thermoacoustic resonators with non-uniform 
medium and boundary conditions. J Vib Acoust - Trans ASME 2011; 133:031012. 

[13] Zhao D, Chew Y. Energy harvesting from a convection-driven Rijke-Zhao 
thermoacoustic engine. J Appl Phys 2012; 112:114507. 

[14] Zhao D. Waste thermal energy harvesting from a convection-driven Rijke- 
Zhao thermo-acoustic-piezo system. Energy Convers Manage 2013;66:87-97. 

[15] Hong B-S, Chou C-Y, Lin T-Y. 2D transfer function modeling of thermoacoustic 
vibration engines with boundary heat-flux control. Asian J Control 2013. 
http://dx.doi.org/10.1002/asic.736 . 

[16] Sujith El, Waldherr GA, Zinn BT. An exact solution for one-dimensional 
acoustic fields in ducts with an axial temperature gradient. J Sound Vib 
1995;184:389-402. 

[17] Hong B-S, Ray A, Yang V. Wide-range robust control of combustion instability. 
Combust Flame 2002;128:242-58. 

[18] Hong B-S, Lin T-Y. Robust observer designs for thermoacoustic dynamics with 
distributed-delay heat input. Asian J Control 2007;9:326-32. 

[19] Matveev KI, Culick FEC. A study of the transition to instability in a Rijke tube 
with axial temperature gradient. J Sound Vib 2003;264:689-706. 

[20] Zhao D. Transient growth of flow disturbances in triggering a Rijke tube 
combustion instability. Combust Flame 2012;159:2126-37. 

[21 ] Culick FEC. A note on Rayleigh’s criterion. Combust Sci Technol 1987; 56:159-66. 

[22] Ferreira MA, Carvalho Jr JA. A simple derivation of the Rayleigh criterion in 
integral form. J Sound Vib 1997;203:889-93. 

[23] Nicoud F, PoinsotT. Thermoacoustic instabilities: should the Rayleigh criterion 
be extended to include entropy changes? Combust Flame 2005;142:153-9. 

[24] Durox D, Schuller T, Noiray N, Birbaud AL, Candel S. Rayleigh criterion and 
acoustic energy balance in unconfined self-sustained oscillating flames. 
Combust Flame 2009;156:106-19. 

[25] Hong B-S. Construction of 2D isomorphism for 2D H CX) -control of Sturm- 
Liouville systems. Asian J Control 2011;12:187-99. 

[26] Hong B-S, Su P-J, Chou C-Y, Hung C-I. Realization of non-Fourier phenomena in 
heat transfer with 2D transfer function. Appl Math Model 2011;35:4031-43. 

[27] Rabenstein R. Transfer function models for multidimensional systems with 
bounded spatial domains. Math Comput Model Dyn Syst 1999;5:259-78. 

[28] Rabenstein R, Trautmann L. Multidimensional transfer function models. IEEE 
Trans Circ Syst I, Fundam Theory Appl 2002;49:852-61. 

[29] Rabenstein R, Trautmann L. Digital sound synthesis of string instruments with 
the functional transformation method. Signal Process 2003;83:1673-88. 

[30] Jin B, Lu X. Numerical identification of a robin coefficient in parabolic 
problems. Math Comput 2012;81:1369-98. 

[31] Xiong XT, Liu XH, Yan YM, Guo HB. A numerical method for identifying heat 
transfer coefficient. Appl Math Model 2010;34:1930-8. 

[32] Hong B-S. A toolkit for developing online operation of gas flow in a confined 
chamber Part 1: model reduction and dynamical measurement. T Aeronaut 
Astronaut Soc China 2003;35:275-87. 

[33] Chu BT. On the energy transfer to small disturbances in fluid flow (Part I). Acta 
Mech 1964;1:215-34. 

[34] Rayleigh L. The explanation of certain acoustical phenomena. Nature 
1878;18:319-21. 




